A Theoretical Study of Heteroj unction and 
Graded Band Gap Type Solar Cells 


Annual Report on 
NASA Grant No. 
NSG 1116 


March 1976 


(NASA--CR-1 46840)- ■ A' THEOEETIGAL STUDY OF N76-21685 

HETEHOJDNCTION AND GRADED BAND GAD TYPE 

SOIAR^ CEDIS Annual Report. (North Carolina 

State UniT,) 63 p fiC $4,50 CSCL 10A Unclas 

. G3/44 25169 


J. R. Hauser and J, E. Sutherland 


North Carolina State University 
Electrical Engineering Department 
Box 5275, Raleigh, NC 27607 




i; ' ‘ 
I 


\ 


APPv i97G 

riEravtrj 




r*'. 
c-‘A 

f' 


V;7f. 



A Theoretical Study of Hetero junction and 
Graded Band Gap T3rpe , Solar Cells 


Annual Report on 
NASA Grant No . 
NSG 1116 


March 1976 


J. R. Hauser and.J. E. Sutherland 


North Carolina State University 
Electrical Engineering Department 
Box .5275, Raleigh, NC 27607 



Table of Contents 


Page 


1. Introduction .• • 1 

♦ • # * . t 

2. ' Theoretical Concepts of Variable Composition Solar Cells 2 

2.1 Features of an Ideal Solar Cell 2- 

2 . 2 Material Parameters Determining Device 

Characteristics 2 

2.3 Limitations of Conventional Cells 3 

2 . 4 Advantages of Variable Composition Materials ........... 4 

3. Device Equations for Computer Analysis 7 

3.1 Equations for Homogeneous Solar Cells 7 

3.2 Modification of Equations to Allow for Variable 

Composition 3 

4. Specific Device Parameter Modeling. .• 16 


4 . 1 Introduction 16 

4.2, Dielectric Constant vs . Composition • • • v"? 

4.3 Band Structure Parameters vs. Composition 17 


4.3.1 Band Gap vs . Composition . 

4.3.2 Effective Masses 

4.3.3 Hole Mobility 

4.3.4 Electron Mobility 

4.3.5 Band Parameters, 0 and 6 

n p 

5 . Results of the Computer Analysis of A1 Ga^ As Solar Cells . . . .• 27 


5 . 1 Introduction 27 

5.2 Material -Parameters for AlASj GaAs, and A1 Ga^ As 27 

5.3 Specific Computer Results 30 


17 

19 

21 

23 

25 


6 , Summary . 58 

List of References -. j 5 59 



1 . Introduction 


This report summarizes research work so far performed on the .theory 

of hetero junction and graded bandgap solar cells.; The 'major objective 

of the work is to investigate the material compositions and, device 

dimensions . needed for high efficiency solar cells. Because, of the 

involved analysis -of -the fundamental equations describing solar cell 

operation, a general ►numerical device analysis program is .being used 

/ 

to study. solar, cells.- 

A major part of the 'initial work has -.been involved , in ’modifying 
an existing silicon. solar cell analysis program .to account for -the 
unique features of graded bandgap, and hetero j unction solar cells. This, 
phase of the work ’has been, completed and the , program is now. being used 
to study. solar cell performance. 

The-, most successful -III-V solar cells have so far been constructed 
in the GaAs and Ga. Al-As material systems i Consequently the work' 
reported herein is concerned with such solar cells . • Th^ most efficient 
solar cell so far evaluated- is an abrupt hetero junction 3ell with a 
pure AlAs layer at the surface with a .GaAs .substrate . • T^e predicted 
efficiency for this cell is slightly, larger than that of a graded, band-r 
gap Ga- A1 As solar cell. 

X ^ 



2. Theoretical ■ Concepts ‘of .Variable •• Composition 

Solar Cells 
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2.1' Features of an -Ideal Solar Cell 

In or.der.to effectively convert . light, energy to electrical energy, 
a solar cell -must .perform several tasks-. First, its 'surface must 
transmit, rather than reflect most of the light striking it. Second, 

■the .crystal lattice of .the cell should absorb the solar photons, 

I 

exciting electrons , out of the valence band and creating electron-hole 
pairs. Third, built-in fields must -separate these excess carriers before 
they can recombine. ■ Finally, the cell should deliver the carriers to an 
external load without, a large 'internal -resist iye loss, 

The efficient -performance of these tasks calls for several basic 
device characteristics, not all of which are simultaneously achievable. 
Outstanding among. these ideal characteristics are: 

1. ' Close optical match -between' cell surface -and external medium 

2. Efficient absorption of -photons in the -solar spectrum 

3.. Carrier generation concentrated in regions of large built-in 
field and 'low recombination rate- , 

4. Large open circuit voltage 

5 ., Small ■ bias current under dark conditions . 

6. Low sheet resistance.. 

2.2 Material Parameters Determining Device Characteristics 

In order to optimize -the solar cell structure, -the material 
parameters that control -the device characteristics must 'be determined .- 
The -first chap?acteristic, a reasonably good optical -match at the • 
semiconductor surface, can usually be obtained with -an -appropriate anti- 
reflection layer that transmits the bulk of the energy at the wavelengths 
dominating the solar .spectrum. 
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The second characteristic, efficient absorption ,of photons, is 
controlled primarily by the value of the ‘band gap, and, absorption coefficient. 
Only those photon energies larger than the band 'gap will create significant* 
numbers of electron-hole pairs, ■ However, when -the photon energy greatly 
exceeds the band gap, much of the excess, will be wasted. ' 

The third characteristic calls, for carrier generation predominantly 
near the. p-n junction -of a cell. Also, carrier generation should. not be 
wasted in regions of high recombination rate, par.tipularly • near the 
surf ace ,, where imperfections cannot be entirely eliminated. Since the 
maxim,van light -intensity always occurs- at the cell -surface, surface 
recombination is -one of the most. difficult problems to solve, especially 
in III-V -solar cells. 

Large open circuit voltage calls for a wide band; -gap semiconductor. 
However,- the wider the band; gap, the smaller the fraction of photons having 
sufficient energy to create electron-hole pairs.- Similarly, a small 
forward dark current depends on a large band! gap with’^he accompanying 
large built-in potential across. the p-n junction. 

Finally, a low sheet resistance can be obtained by using a thick 
surface layer above the -p-n junction.- However, the thicker the. surface 
layer, the greater the attenuation of light reaching, the r junction 
depletion region, the most desirable area of optical generation, 

2,3 Limitations of Conventional -Cells 
The fact that conventional solar cells are usually constructed .of ' a 
single material (usually silicon) severely limits the degree to which 
the -cell characteristics can be influenced. One. of the most, crucial 
material .parameters, the band gap, becomes- fixed as soon as the cell 
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material is selected. A constant band gap results in a monotonically 
decreasing optical generation rate. In other' -words , a homogeneous. solar cell 
concentrates carrier generation near the surface,, where ’the recombination 
rate is highest, and the -built-in field from the -p-n junc;tion is relatiyely 
weak. If the p-n junction is moved very near, the .surface to separate 
the .electron-hole pairs more efficiently, the surface layer sheet 'resistance 
becomes unacceptably high. 

In addition,- silicon has a 'band gap. (abput 1.1 eV) which is ^ too small , 

to efficiently use -the photon energies that dominate, the ;Solar spectrum 

[l-,2]. Theoretical calculations have estimated that a .band, gap value .of 

/ 

about 1,4 eV .should 'lead to , the most efficient carrier generation by 
sunlight. Too few solar photons 'have sufficient energy -to -bridge ' larger 
gaps . ■ 

In addition to more efficient ' carrier generation, .a larger band: gap. 
semiconductor leads to a higher, open circuit -voltage [2,19], and lower 
dark' current , both "mentioned earlier as ..desirable characteristics. 

In short, conventional -homogeneous silicon ■ solar cells leave much to be 
desired in .theoretical efficiency. Their main attractiveness. stems -from 
the advanced state of silicon technology. 


2.4 Advantages- of Variable Composition Materials 
Clearly., ,the variation of material composition with depth, available, 
in III-V semiconductors should > permit -improved -control over- the .distri- 
bution of carrier generation. The band diagram in Fig. ..2.1 shows a wide band - 
gap (E surface layer covering a smaller band, gap (E ) .substrate. 

If Eg^ is chosen large .enough, the -surface layer -can serve ‘as .a window, 
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transmitting most photons to the p-n, junction where the band gap is , 
smaller. Thus =, carrier generation can, be ‘concentrated in, a region of 
large built-in field, .far. from surface recombination sxates. If the 
concentration of interface states in the regions of variable composition 
can be kept at acceptable levels, then minority carrier collection should 
be accomplished with much greater efficiency than - in homogeneous solar 
cells o 

In addition, a .thick surface layer can.be used to reduce sheet 
resistance, without • severely attenuating the light intensity reaching 
the- p-n junction. Also-, band gap grading can be used to create 
additional built-in fields that aid -carrier collection. [ 3,4] , In fixed 
band gap, devices, built-in fields must 'be .created- by doping gradients. 
However, if the band gap can. also vary,' independent, effectiye built-in 
fields can be. made to operate on electrons and on holes; 


E 


no 


1 

q dx 


po q . dx 

-Finally, the' cell shown, in -Fig., '2.1 should have a .larger- open -circuit- 
voltage than a comparable homogeneous ' cell 'made of .a material with the • 
smaller band gap, same time, such a variable gap cell makes • 

effective use of -the lower -energy photons that, would be .lost. -to a 
homogeneous cell 'made of a material with- the larger band gap , . 

Thus, the variable composition solar cell appears ‘to be a step 
closer to the ideal cell described in Sec,. 2.1. • However, the, complexity 
of the trade-offs involved defies .-simple .approximations .. For this • 
reason, a computer analysis program is desirable to give a better estimate 
of the .performance improvement 'that -can be expected. 
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3 . - Device Equations for Computer .Analysis • 

3 , 1 Equations for Homogeneous Solar Cells 

In order to adequately -analyze variable band gap solar cells, 
several device phenomena .must 'be incorporated into the- mathematical 
equations used to, model the .cells. Many effects included in the present 
work are important in ordinary homojmction solar cells, - Among these 
are; • 

1. Drift and 'diffusion currents 

2. Position dependent doping 

3. Doping dependent mobility 

4. ‘ Optical carrier generation 

5. Bulk generation-recombination effects - 

6. Surface recombination effects 

The '.introduction of spatially varying composition demands that • 
additional factors be accounted -for: 

1. ’ Position dependent band gap 

2. Position dependent electron affinity, 

'3, Built ;in fields ,due to. a varying band 'gap 

4-, Composition dependent refractive index at the surface 

5.' Other position dependent material parameters such 'as: 

a, mobility 

b , dielectric constant 

c, optical .absorption coefficient 

The -'basic semiconductor equations used for homoj unction devices 

t * 

can.be modified to model variable composition device behavior. The 
unmodified equations in one dimension are '[5-9]: 
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Transport : 


J = 
n . 


« „ 9n . 

qy nE -+ q -D i 
^ n ^ ■ n 8x 


(3a) 


J 

P 


= qPpPE - q Dp 



(3.2) 


Continuity;, 


8n _ 
3t 

, 3J 

- U + + G , 

q 3x e 

(3.3) 

3t 

1 

- U - — . ^ + G , 

q 3x . 

(3. -4) 


Poisson's -equation; - 




Auxilli.ary equations: 


U = 


np-n^ 


T (p+p )+T / (n+n-T* 
no ^ po 1 

(Shock ley-Read-Hall model, with single trapping level) 

Epn~E^ 

n = exp[ (non- degenerate) , 

p = N exp[' i"v^^J^ ~ ] (non-degenerate) , 

V Ki 


N = (net doping) 


(3.,5) 


(3.6) 


(3.7) 

(3.8) 

(3.9) 


„ (N>0 for n-type). 

This system of equations is valid only for a device made of a 
single semiconductor material doped with enough impurities to create 
the desired doping profile. The -electron energy band structure and. its 
associated, characteristics must be the same throughout the device. 
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The .existing computer device analysis program for homogeneous solar 
cells, [9] depends on the fact that ‘the preceding group of equations . can 


be reduced, to three equations in the three unknowns'.’ and <1)^, where, 

i|/ is the electrostatic potential and'ij)^ and are the Fermi potentials, 
defined by: 


^n ■ q Fn’ 


(3.10) 


= -,-E„ . 
p q Fp 


(3.11) 


It can be shown that .the same basic approach using the same three 
variables, can be used- to analyze variable, composition solar cells. 

The modifications required to account .for material variations are. 
discussed in the next section. 


3,2 Modification of Equations , to 
Allow for Variable Composition 

In order to rewrite the device equations so that they apply to 
variable ‘composition devices, it is useful' to recall the basic principles 
from which they were derived. In .particular, Equations . (3 .1) and (3.2). 
can be written as: 


J 

n 


d E 




Fn 


dx 


(3.12), 


. d E 


J • = y p 
P ■ P 


12 . 


dx 


(3.13) 


For a homogeneous material, it can be seen that Equations (3.1) and ‘(3. 2) 
result from Equations ’ (3,12) and (3.13) in the following manner. ■ From 
Figure 3,'l;, 


^c “ ‘^'^o ■ ^ " X-c» ■ (3.1i^) 

and . 


E 

V 


- 

o 


q 'I' - 


g 


(3.15) 
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From equations (3.7) and (3,8), 


E„ = E +kTAn(^), 

Fn c N 


(3,16) 


^ 'V 


(3.17)- 


Using Equations (3.14)- and (3.16), Equation (3,12)- can be rewritten as:. 




(3.-18) 


rp . kT dn kT ^c-, . 

^n^ 0. dx • n .dx ’ N . dx 

c 

And, using Equations (3.15)' and- (3.17), Equation (3.13) can be rewritten 


as: 


= P„P It 

V' 


■p dx “"""o “c ~g 

dE, 


(3.19) 


dN. 


^pP dx 'dx p dx N dx - 


For a homogeneous material, x ) N , and N are constant throughout 

O g O V 

the device.. Therefore, Equations (3,18) and (3.19) simplify to Equations 
(3.1) and (3.2).- But for a variable composition ‘structure, Equations 
(3.1) and, (3. 2) must be replaced by Equations (3.18) and (3,19), or by the 
simpler forms of Equations (3.12) and. (3'. 13). 

The continuity. Equations , (3.3) and (3.4), apply equally, well to 
heterogeneous , or homogeneous materials. Poisson's -equation, however, 
fails to reduce to. the simple expression of equation (3.5) when -composition 
(and theriefore .dielectric .constant,, e) is allowed- to vary with position. 

A .more. general form can be ■> derived as follows. 

Basically Poisson's equation in., one dimension is 
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<i(e E) , 

-*r“ = 




= _a.( 

,2 e ^ ' dx ax 


(3. -20) 


Thus,- when the dielectric- constant becomes a. function of -position, 
Equation (3*20) must replace Equation (3.5).. 

The carrier concentration equations demand no modification,- but 'they 
can be written -in . a more convenient form for the purposes of computer 
analysis, Using Equations (3.10), (3,11),.- (3.14) and (3.15>, Equations 


(3.7) and (3.8) become: 


. . kT. „ 


n = n.^.exE [J. 5 ^)] 


1’ = "12 .""J" ^ - "i * -'•'o 

' « r 

where n^^ ~ intr;lnsic carrier-- concentration of material ,#2 in -the alloy ., 

The 'composition dependent , (and. therefore position dependent) parameters 


can be -grouped into two .terras 


^c kT ^c 

0 ■ ^ 

n q ""o q n^g 


(.3.21), 


0 -= _ 1 +E ) + .,;; + ]Sl' to 

p q ^c =g ^o q. n ^2 


(3.22) 


Then, the, carrier density equations can be written concisely as 


n ■= n^2 ^n-^ » 


(3.23) 


p = n ^2 


(3.24) 
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These expressions are'analgous to the forms used in the, computer 
analysis of homogeneous solar cells [9]: . 

n = n^ exp 1 

p = n^ exp 

The parameters j 6^ and 9^- depend only on composition, and must 'be r 
externally supplied as, functions of position after , the composition 
profile is specified. The. use of thes.e parameters , is dis.cussed further 
in -Sec, 4,3.5, t 

The' other amcilliary equations remain the same whether or not ■ 
composition is spatially .varying. However,, it is .-important to note 
that several 'Variables no longer can be assumed- to be constant. In -fact, 
the material parameters . , y^, e, n^, n^ and p^ must be given as known 
functions of composition (or position, when thsi composition profile is 


Thus,, the .complete one dimensional mathematical 'model for -variable 
composition solar cells consists of the following system of equations: 


dE_ 

■J = y n.'— 5 — , 
n n dx ■ 


(3.25) 


J ' 

j . = • y ' p 

p p^ dx ’ 


(3.26) 


is. " G U + ^ ■ 
_ _ - U + - 


(3.27) 


^•=6 - U - i ^ 

8t e q dx’ 


(3.28) 


= _ a (p-ntN)- 


(3.29) 
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n 


np-n. 


U = 


T (p+p, )+ T (n+n.' ) ’ 
no ^X' po 


ti.o e^P (ip-<j>„+e„)] 


12 


n. n 


P = n.2 exp [^.(^,^-*+6^)], . 


with, auxilliary equations: 


Tn 


q. 4' 9 

■a Yjj9 


E_ = - q <t> , 

Fp ^ '^p’ 


and 'with 'the, following parameters .given as functions, of position: 
N = N(x), = net_ doping profile 


C = C(x), composition profile 

(for ex^ple,' C .= mole fcaotion, x, 
of -AlAs in an A1 Ga., As. solar cell) 

= y^(x), electron .mobility 

Pp = Pp.(x), hole mobility' 

e = e(x), dielectric constant 

n. =,n. (x), intrinsic carrier concentration 
11- 

= n^(x)-, trapping center parameter for electrons 

Pi ='p^(x), trapping center parameter for holes 

t = -T (x), excess carrier lifetime in p"*" material 

T ='T (x), excess carrier lifetime in n"*" material 
po po ’ • • 

E = E (x), band- gap 

o ^ a 

0 = 0 (x) , conduction band parameter 

0 =6 (x) , valence band parameter 

p. p . r . 

G ='G (x), optical generation rate 
$ 0 


(3.30) 

(3.31) . 

(3.32) 

(3.-, 33) 
(3.34) 

(3. .35) 
(3.36) 

(3. -37) 

(3.38) 

(3.39) 

(3.40) 

(3.41) 

(3.42) 

(3.43) 

(3.44) 

(3.45> 

(3.46) 

(3.47) 
(3.48> 
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Magnetic effects and thermal ■, -gradients are assumed to be negligible. 

The modeling of several -of the above parameters -will be discussed in 
the following sections. 

This t apparently unwieldy mathematical ■ system can actually be reduced, 
to three equations , in three- unknowns in exactly , the same manner used 
for homogeneous solar cells in previous computer models, The fact that 
composition (and -the associated, material parameters) is allowed to vary 
through the cell simply adds some complexity, to the detailed computations 
without altering the method -of solution. For this reason, it has been 
possible to modify an existing computer analysis program designed for 
silicon horaojunctipn sol^ cells f9], in order to permit analysis of 
variable composition cells. A major share of the work up .to -the present 
time has been devote.d to thdse detailed modifications: However, the ■ 

modified program is now essentially complete and .-has been used to examine 
several hetero junction and, graded 'band gap structures. The results of 
these studies will.be discussed in "a later .section. First, however, the 
modeling of several of the material parameters deserves some, explanation. 
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4.- Specific Device Papameter Modeling 
4,1 Introduction 

The 'modeling -modifications discussed in the following sections have, 
been. used -to convert a computer program designed for -a single 
homogeneous 'material (silicon) into a program capable of. examining 
a spatially varying alloy of two compatible semiconductors . Given the. 
material properties of the ,two semiconductors , the^ program . can theoretically 
analyze the 'performance, of a solar cell of any specified composition profile. 
For example, solar -cells made of the .alloy, A1 Ga, As, have been analyzed, by 
providing the material properties of AlAs and GaAs, along, with the desired 
composition profile, i;e.. the ;mole. fraction ;Of AlAs, x, vs. -position.. 

Because .experimental data concerning the, properties of semiconductors and 
their alloys is limited, reasonable approximations must be used to predict 
the properties of arbitrary alloys.- Whenever possible, these approximations 
have been. checked against experimental results and, adjusted to improve the 
agreement with available data. The most important approximations are 
described -in the following .sections . 



(4.1) 
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Exactly the' same form is used to estimate the high frequency relative 
dielectric constant, for an allpy: 


- ~1 - e, 1 

i+-2[ci:--^) + ( i-c) : 


= 






'h2 


1-C(r^)-(1-C)(-^) 




£ +2 
®h2 


(4.2) 


4,3 Band Structure ■ Parameters ■ vs . Composition 

The original .homogeneous solar cell program was based on the -assump- 
tion that the material being studied had a band structure with a. single 
dominant ' conduction band,-valley , either direct or indirect. However, it 
is possible -for a variable composition solar cell to be ‘ constructed of an 
alloy of a direct band gap material and an indirect band gap material 
ClOjll], Depending on the relative conduction band -minima of the two 
materials and. the alloy composition, . it is .possible .that -two conduction 
band valleys . (one -direct and one, indirect), significantly affect carrier 
behavior . 

In-order to allow for such a situation, the- modified program permits 
the .specification of two valleys, each with its own . composition dependent, 
band gap and 'its own characteristic effective mass and mobility. The 
following sections describe the techniques used to determine the -resultant 
band structure .parameters (such as electron and hole .mobilities), of an 
alloy, when the band. structure parameters for each of the two components 
are given. 

4. -3.1 Band Gap vs. Composition 

Referring ,to Figure 4.1, the resultant direct and indirect band gaps are 
estimated empirically from the gaps of the component materials. 




Figure 4.1 Electron energy, E, vs. wave vector, k, for a) indirect 
semiconductor, material 1, b) direct semiconductor, 
material 2,. c) alloy of 'materials 1 and 2 
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Experimental studies have shown -that', when two • materials having • direct 
gaps, ^gd2’ alloyed, th& resulting direct gap has .-a quadratic 

dependence . on composition: 


E , = a + bC + E (4,3) - 
gd gd? 

where C ,= mole fraction of material 1, and a and b.are adjusted f or - 
closest fit; to the. experimental data. For A1 Ga, As,, the indirect 
gap was, estimated by. linear interpolation: 


E . = E •.. + (E .„-E .,)(1-C) 

gi- gil gx2 gil 


(4.4) 


The resultant alloy band gap is,, of course., given by: ' 


E = Minimum of(E,,E.) (4.5) 

g ’ gd. gi 


4,3,2 Effective Masses 

The next task is, to determine :the effective .masses for holes and 
electrons. ‘Since the valence -band ‘in. each of the •two , semi conductors 
is .assumed' to have a single dominant yalley with a characteristic effective 
mass, it is simply, necessary to .interpolate' between the ^ two 'effective 
masses to find the resultan't; hole effective mass 'for the alloy. The 
following interpolation form was used- [20]: 


Jv = ^ + 

m * m*i m*„ 
p. pi. p2 

where m*, = hole effective mass in material 1 

pi,. 

ra*_ = hole ■ effective mass in material 2 ' 

p2 

m* = hole effective mass in alloy 
P 

C = .mole fraction of 'material 1 


(4.6) 
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The same form was used to find the electron effectiye mass for the 
direct valley: 


C ^ 1-C 
-j : — + 


cd cdl cd2 


(4,7) 


where m*^, = conduction band effective mass in direct valley -for material 1. 
cdl 

m*^2 = conduction band effective mass in direct .valley for material 2; 
”*cd ” conduction' band ' effective mass- in direct valley for alloy 
and likewise for the indirect 'valley: 


m- 


m»’ 


1-C 

m* . ^ 
Cl 2 


(4.S) . 


ci cil 

Now, 'it will prove useful to, define, if .possible, a. single effective 
mass for all electrons. This can be .done in the. following way: 


n = exp[ 


^cd 


kT + ^ci 


Ep -E . 

Fn ci 


kT 


•] • 


where = direct valley effective density of states 

=: indirect valley effectiye density of states 
E^^ = . direct valley minimum 
= 'indirect 'Valley mln-lmum 

It is desired -to obtain the following form: 

n-= expC-^2g^] . 

where = effective -density of states -for entire conduction band 

E = -absolute minimimi of conduction band (the lesser of E', and 
c ... ^ 
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Since N , = K 

cd cd 


N . =,K 

Cl Cl 


3/2 


■kT ir a3/2 

N = K m-* • 
c n 

where K is a proportionality constant 'and m* = effective mass for all 
electrons . 

Then, 

_.i3/2 r^Fn'^o-, 3/2 

■»". exp [— jgr-] = exp C 

, 3/2 r^Fn ^CIt 

+ mj. exp C— ig— ] 


1 r X 3/2 , 

11^ = exp( ^ ^ 


E -E E -E • 

c '"cd,, ... 3/2 , c ci,-,2/3 

)+m-. exp(— — )] '' 


or, in -terms of band gaps : 


m' 


n = exp(- \ - /-0 -h exp( -^ -- ) ] 


n 


kT 


Cl 


kT 


(4.9) - 


where E ,, E E , m * . and m’*. are given by Equations (4.3),- (4.4), (4.5), 

gQ 5 CO ‘ OX 

(4.7) and (4.8). 


4.3.3 Hole Mobility 

Thus -far, the effective masses for holes and electrons have been 

determined as functions of composition when the. material parameters of the 

two alloy components -are given. Since the device equations used to model 

the -solar cell require mobility data, instead, of effective mass, the final 

step is to determine hole and electron -mobilities . 

! 

Mobility depends not only upon host material composition, but upon 
the doping level as well. In order to calculate hole mobility for 
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arbitrary doping and. composition, an -.empirical formula, ^p 2 »- ' 

constructed to approximate the doping dependence of .hole mobility for 
one of the alloy components, to be referred -to , as material 2: 


where N = total impurity concentration 
Now, because 


(4. .10) 


q<T > 
P 

V = — 
p m- 


<"p2^ 


<T > 

P 

m* 


(4.11) 


it is apparent that if <t^> and m* can be: specified as functions. of 
composition, then a reasonable approximation to hole . mobility for any • 
doping level and composition can be -made. If mobility is controlled 
primarily by polar optical phonon scattering, then ”^Tp^ can 'be -written 
as [2l3; 


K 


p r~sr , 1 1, 


/m^ (t ) 


(4.12) 


p e 


whei;e K is a proportionality constant for all materials . 

= high frequency relative :dlelectric -constant ' 

= low frequency relative dielectric constant 

Since .Equations (4.1),- (4,2) and (4.6) give e^, and m* as .* 
composition varies, the modeling of hole mobility, is virtually complete. 
Combining Equations (4.10), (4.11) and (4.12): 


Vp(N,C). 


■e 3/2 1 1 

._ ^p2 ^^h2.~ £^2 


) • 


,.3/2 V 1 , ■ 1 ^ 

m-' ( ) 

P . £/ 


(4.13) 


'h- 


where ni* 2 s ^£ 2 ^ ^h2 hnpwn parameters of material 2. 



23 


fp2(N) is an empirical ‘function relating "hole mobility and doping 
in material ■ 2, 

and' and m* are given as .functions of composition by. Equations 

(4,1), (4,2) and (4.6), 

This , completes the ;specif±cation of hole mobility as a function of, 
doping, N, and composition, C. 


4.3.4 Electron Mobility 

The, ihodelling of electron mobility is more . complex than hole mobili-ty 
if an alloy of -an .indirect .gap semiconductor .(material 1) and 'a direct 
gap semiconductor (material 2) is considered. The .following method 
depends -on the assumption that the ;dopiiig dependence of mobility for the 
for the Indirect valley of the alloy can be inferred firom data on material 1, 
and' the dependence for the direct band of the •, alloy can be -inferred - from - 
data on material 2. Treating the direct and indirect 'valleys Independently, 
and using exactly the. same technique as. described for hole .mobility., a 
direct valley electron mobility, and an indirect valley electron 
mobility, can.be -determined as: 


y^(N,C) . 


y^(N,.C) 


-p \ A 3/2 [ 1, 
”Sd2. 

. J.) 

CCl C i> 

h K. 


nl ci2 ■ 




m*. 

Cl 


3 / 2 ^^.^ 



(4.14) 


(4.15) 
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wher;e 


_ - A 

HI". 


cil’ 12* IV 


^h2 ^hl knovm parameters of material 1 


and 2, 

is an -empirical function relating electron mobility, to 
doping for material 2 

f^^(N) is an empirical function relating electron mobility to doping 
for material 1 

and Ej^, m*^ and.m*^ are given -by Equations (4.1), (4.2) ,■> (4.7) and 

(4.8). 


Now, the final task is to merge ‘these- two mobilities -into, a single 
resultant electron .mobili-ty. Such an ’'effective" mobility can be found, 
by weighing the direct aqd indirect mobilities -by their respective 
electron populations. These can be determined by the following procedure: 
define = fraction of -free electrons in direct valley 
n^ = electron density in direct valley 
n^ = electron density in indirect valley 

Then, 


R , = ---- 
d n' +n. 
d - 1 


„ ,^Fn 


E -E E“ ~E 

, Fn cd, , „ , Fn ci. 


kT 


kT 


or, - 


^d- = 


m“. 

l+(-^) 

“cd 


3/2 


E , - E . 


(4. ,16) 
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where the composition dependent .parameters are given by Equations (4.3),- 
(4.4), (4.7) and (4-,8).- Finally, the resultant electron mobility is 

4^'(N,C) = (4.17) 

where Equations (4.14) , (4.15> and (4.16) give -p^, and in terms 
of composition. 


4.3.5 Band Parameters, 6 and 0 

L_2 E. 

In order to model the composition dependence of the band parameters, 

6^ and 6^, some knowledge of the variation of -electron affinity, Xj,s 

band gap, E , must 'be available. It has been found to be .advantageous to 
S 

select the potential reference , ip , so that 6 =6 =0 when the mole 

o n p 

fraction of material 1 is zero. In other words, 6 and 9 are zero in 

n p 

regions "of the solar cell ..consisting entirely of material 2. Then, 
from Equations (3.21) ^nd'(3.22): 


o q q n^2 


(4.18) 


Y + E N 

♦ = (^) 


"i2 


(4.19) 


The fact that these . two conditions are equivalent can be seen by 
manipulating the well-known expression for intrins;ic carrier concentration: 


2 - 


E -E 

F ^c. 


E -E^ 

V F, 




v2 


kT 


kT .en 


'g2-' 


N .N _ 
c2 v2 


(4.20) 
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Novj, inserting Equation (4,18) into (3.21), 
gives; 


e 

n 


Z -X „ , _ N 

( 0 ) 

q. q »^2 


and (4,19) into (3.22) 


(4.21) 


0 

P 


^c2-^c 


P —H N 

q q . 


(4.22) 


Finally, in terms of effective mass: 


e 

n 



]<1 

q 


m-* 


In (-^) 
111'* 
n2 


(4.23) 


0 

P 


-(Ax^+AE^) 

_ 


■.+ 


3 

2 


kT 

q 


m* 

■tn (— § — ) 

™p2 


where Ax^, = X^ - X^2 



= E -E 


g 


g2 


(4.24) 


Equations (4.23) and (4.24) are the forms .used to compute the band 

parameters for an arbitrary alloy of two materials. The band gap and the 

effective masses are given as functions of composition by- Equations (4.5), 

(4.6) and (4.9). The modeling of electron affinity vs. composition depends. 

on the experimental data available for the alloy being examined. A t 3 rpical 

case of hZ Ga„ As is. discussed in Sec. 5.2. 

X 1-x 
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5. Results of the Computer Analysis of 

A1 Ga. As ■ Solar .Cells 
X . 1-x 

5 .1 Introduction 

At th^ present ■ time the program modifications .are essentially 
complete and the program has-been used to. examine several hetero junction , 
and variable band gap solar cells made of A1 Ga. As • As indicated .in 
previous -sections,' the program requires considerable, input .-data to 
define material "parameters , The ..following section .lists the .material , 
specifications and assumptions , upon .which the .results for A1 .Ga. As 
are based, 

5.2 Material Parameters for AlAs, GaAs, and A1 Ga. As 

JC «l«*rX 

Material ^^1:. AlAs 
Material #2; GaAs 
Dielectric constants [22]: 

low frequency; =-10.9 

^L2 " 

high frequency: , =..8,5 

Effective Masse.s C29];" 

• holes; m*, = 0.85 -m 
pi ' o 

m*;, = 0.-6S m 
p2 , o 

electrons ; . 

direct valley; m*,. = 0,126 m 
• cdl ' o 

rn’K- = 0,0636 m 
cd2. ■ o 
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indirect valley: , in*.^ = 0,37 m 
^ cil o 

m*:_ = 0.39 
ci2 o ■ 

7 —3 

Intrinsic carrier concentration for-GaAs: n.„ = 1,1x10 cm 

i2 . 

Temperature =.300°K 

Trapping Center: Single trapping center assumed to be at center of band 

gap regardless of -alloy composition: 


n, = ,p. = n. 
11 1 


Mobility vs. doping, N: 


Holes,- data for GaAs -used [12]:- 


= f ^(N) =• 


360 cm /v.s 


Clt(3.17xl0-"^cm^)N3^'263 


where 'N = impiirity sites/cm , 

Electrons, indirect valley (AlAs data used) [23]; 


^ (•vj\ _ 165 cm /v*s 

^nl " nl^ "" rTTTTTTrlT 3.,,,0.13 
Cl+(8,lxl0 cm )N] 

Electrons, direct valley (GaAs data used) [12]: 


y = f (N) = Z20.Q._.9”L /lls.: 

Band Gap vs. composition: 

AlAs: direct gap = =.2,95 eV 


indirect gap = 2.16 eV 

overall gap = E , = - E . ^ 

^ gl gil, 


GaAs: direct gap = ^^^2 ~ 1*^^® 

indirect gap = - 

overall gap = E^^ = 


[ 17 ] 
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Alloy direct gap, useli empirical formula [17]: 

E , = 0.468C^ + 1.042C + 1.^39 
gd- 

( where C -s mole fraction of AlAs) 

Alloy indirect gap* used linear Interpolation: 

E . = 2.16 + 0.29 (C-1) 

gx V 

Electron affinity vs. composition [24]:' 

(Only 'Ax^ = “ Xq 2 required for the analysis. 

X ^2 assigned a value. See sec. 4.3.5) 

Lifetimes -[27]: t = 5.3 x 10 ^ sec 

no 

-Q 

T = ■ 8. 5 X 10 sec 
po 

(a first ‘Order model to be improved later) 

Surface recombination velocity: 

5 ■ 

S = 10 cm/seq affront surface 
S = 0 cm/sec at back surface 

2 

Solar power density at solar cell surface: 135,3 mw/cm , [25]. 

The-, optical generation rate, G (x) deserves some explanation. 

t 

G^(x) was calculated by an independent program similar to the one used 
for homogeneous . solar cell analysis ^9],. This modified- program 
accepts an arbitrary composition profile of an alloy of -two materials 
covered by an antireflection layer of SiO, of specified thickness.. 

The -program can be -run in one. of -two modes: standard mode, and efficiency 

mode. By using the .efficiency -mode, it was, found that an antireflection 

o 

layer thickness of 700 A provided the closest match to the solar- spectrum 
(AMO conditions) for any mole fraction of AlAs at the -surface. . The 
program has been designed to permit specification of -either AMO or AMI 
irradiance conditions. However, at the present -time, only AMO conditions 


have been examined. 
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The experimental -data on the refractive index of GaAs, shown 
in Figijre -5.1 [12], -was used to determine the transmission, absorption and 
reflection coefficients at the surface. For lack of better experimental 
data, the. refractive index ‘of ,A1 -Ga. As was assumed, to be. the same 
as that for GaAs, This assumption can be easily Improved as .soon. as 
better data becomes available. The absorption coefficient for A1 Ga. As 
vjas calculated by interpolating between data for GaAs [16], and AlAs [13]. 
The results of this interpolation are shown in Figure 5.2 for ‘six values 
of the mole fraction of AlAs. 

Using the same theoretical techniques as in the homogeneous solar 
cell analysis [9,], was calculated for each composition profile 

and stored for use in the main solar cell analysis- program. Typical 
results are .shown in the next section. 

5.3 Specific Computer Results 

Figure 5.3 illustrates- the structure and.dopipg profile common 
to. all the A1 Ga., As polar cells examined so far . Although only 
abrupt junction n-on-p devices .have been examiped,^ the .analysis 
program permits specification of a Gaussian doping profile, and a. 
p-on-n configuration, qs -well. 

Since -the analysis program was completed, seven . significant 
composition profiles have been examined. These .profiles are listed 
in Figure 5. -4 alongside the cell numbers assigned to each one. It 
should, be noted -that -although the p-n junction of each , cell .coincides 
with the termination .of composition grading, ther.e is .no requirement 
that this be so. The possible advantages of placing the- p-n junction 
Ipside .or completely .outside the -region of composition grading will 
be studied in futiire runs, ■ 
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Figure 5,3 Structure 
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Figure 5,4 Device composition profiles 
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In order to analyze a solar cell,, the .computer .program requires the 
specification of device structure . along with the desired values, of 
terminal voltage; The -results are obtained -in both tabular and 
plotted form. For each value- of terminal 'voltage , the corresponding 
current density, power density and efficiency are calculated along with, 
the values of electrostatic potential, electron and hole current ..densities , 
and electron and hole concentrations at each of -approximately 1000 position 
values,. ' Figures, 5.5 through 5.13 present the plo-tted data for cell 2 , 
as -typical of that which has been obtained for each of -the devices ‘that 
ha? been. examined. In .the plots covering the- entire device (125 yM), 
the- solar cell surface is on .the right-hand s.ide, while in the -expanded, 
plots, covering only the top layer (1.5 yM) , and in the plot of the 
generation rate,, the , cell , surf ace is -on the lefthand side. It -should, 
also be noted that Figures -5.6 throi:igh 5.9 plot the absolute ■magnitude 
of the -current densities, so that sign reversals .appear -as -spikes .- 

The -effects of composition grading are clearly evident 'in .several 
plots. Figure 5. .5 shows the built-in potential resulting from , band. gap 
grading as .mentioned in Section 2,4. The slope of, these curves gives .the 
built-in field. The -.effect 'of band. gap grading on .the optical ' generation 
rate. .can be seen in Figure 5.13,. where the -reduced slope in .the- surface' 
layer 'is due to the absorption of progressively smaller values of photon 
energy as the band gap shrinks .with -depth, 

The- transition from an, indirect to a direct 'band 'gap structure 'has 
a marked .effect on device characteristics. Since, for ,A1 *Ga, As, 
this transition occurs at about 0.4 mole fraction of AlAs [18], cell 2 
consists entirely of a direct 'gap material, and no transition occurs. 
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Figure 5,10 Electron and hole concentrations across cell 2 







Cell. 3, however* has a linearly graded. region from 100% AlAs at .the 
surface to 0% AlAs at a depth of 1 uM,. so that the 'direct-indirect 
transition occurs at -about 0.6 pM below the surface, , Figure 5,14- 
shows the pronounced . change in the built-in potential ,at the point of 
transition* while the composition grading itself is linear. The 
corresponding energy band diagram and carrier concentrations for cell 3 
are shown in Figiires 5,15 -through 5.17. 

The generation rate profile of cell 6, shown in Figure 5.18 
graphically illustrates the improved control over carrier generation 
rate , distribution made possible by composition variation (see Section 2.4). 
This device consists of an AlAs .layer on GaAs. The surface layer acts as 
a window to separate the high recombination rate at the surface from the high 
generation rate region below the hetero junction. The improved collection 
efficiency is responsible for making cell 6 the most efficient -solar cell 
structure of. the seven examined so far. 

The relative performance -of five of the most significant cells is 
depicted in Figures 5.19 -fhrough 5.22 with.no correction for ohmic 
contact area. Although cell 2 has the .largest short circuit current,' 
cell 6 has the largest open circuit voltage and the highest peak 
efficiency, The main advantage of cell 6 over cell 2 is the greater 
concentration of carrier generation within the depletion region 
(compare Figures 5.13 and 5.18). The built-in field, due to band gap 
grading in cell 2 fails to collect carriers efficiently enough to 
compensate for the fact that the maximum carrier generation rate occurs . 
at the surface. Thus, cell 6 has a peak efficiency of 19.19% while 
cell 2 has >a peak efficiency of only 17.92%. 
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Figure 5.14 Electrostatic potential near the surface of cell 3 
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Figure 5.17 Electron and hole concentrations near the surface of 
cell 3 
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Figure 5.20 Current density vs. voltage characteristics for 
cell 1, 2 and 3 
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Figure 5.22 Current density vs. voltage characteristics for 
cells 2, 6 and 7 
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Although the performance -data referred to so far is useful for 
the comparison of cells that have been analyzed by the computer, 
some correction for contact area must be made before a realistic , 
estimate of practical efficiency can be made. Table .5,1 summarizes 
the . significant performance parameters of the 7 cells analyzed up to 
the present time. The only correction for contact area appears in'the_ 
last column 'of , Table 5.1, where the idealized efficiency .has been 
multiplied by a factor of 0.87 to correct for an assumed 13% metal 
contact coverage. Naturally, such a first order correction fails to 
allow for the fact 'that a given , contact geometry does not alter the ■ 

J vs. V characteristics of different devices .in the, same way. Also 
this simple correction factor does not 'account for the sheet resistance 
of the surface layer with a finger contact geometry. A more accurate 
correction for .contact area and series sheet resistance can be 
expected in future analyses. 

Figiire '5. 23 compares the. results of an .analysis by-Hutchby [26] 
of several devices having linear composition grading , with -the analysis 
of similar devices (cells 1, 2, and .3) by .the variable composition solar 
cell analysis program. The device par ameters • of cells 1, 2,. and 3 closely 
match those of the upper curve (see Sec. 5; 2) in Hutchby's data and it 
is apparent that while -the computer predicts- slightly lower .efficiency 
at each value of composition,, the general dependence of efficiency on 
surface composition appears to be the. same in both analyses . 







































Efficiency (Percent) 




Figure 5.23 Efficiency vs. surfac| mole fraction of AlAs, at the surface. .O : S= 1x10 cm/sec. 

L q= 2.1 ym; □:S=lxlO cm/sec, L =2.1ym; <>:S=lxlo5 cm/sec, L =0.52 ym; 

A:S=lxlO^ cm/sec, Lp^=0.52 ym [22]. Xrcells 1,2 and 3 with characteristics '' 
comparable to O . 
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6 . ! Summary 

As the preceding results indicate,- the 'variable composition solar 
cell analysis program predicts efficiencies that -are in agreement with ■ 
qualitative estimates ,of solar cell performance. The 'shift of peak 
optical , carrier generation away- from the surface to the [junction 
depletion region, made possible by position-dependent composition, has 
led to the ejected increase in efficiency oyer that of homogeneous, 
solar cells. 

Future work will include: 1) an examination of ’the possible 

advantages of placing .the- p-n junction above or below the termination of 
composition grading, 2)- a search 'for an optimum composition profile and 
the corresponding .doping profile, 3) improved modeling of the .dependence 
of lifetime -on position, and 4) an examination of solar cells- using alloys 
other than M Ga. As. It ds hoped that -the analysis, of various materials 
and structures .will -lead. to a better understanding of the processes 
controlling variable composition solar cell performance. 
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